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Abstract 

We investigate the behavior of systems of interacting diffusion processes, known as 
volatility-stabilized market models in the mathematical finance literature, when the number 
of diffusions tends to infinity. We show that, after an appropriate rescaling of the time 
parameter, the empirical measure of the system converges to the solution of a degenerate 
parabolic partial differential equation. A stochastic representation of the latter in terms of 
one-dimensional distributions of a time-changed squared Bessel process allows us to give 
an explicit description of the limit. 



1 Introduction 

Recently, Fernholz and Karatzas Q have introduced two types of systems of interacting diffu- 
sion processes, the volatility-stabilized market models and the rank-based market models, in the 
context of stochastic portfolio theory. Both of them serve as models for the evolution of capi- 
talizations in large financial markets and incorporate the fact that stocks of firms with smaller 
market capitalization tend to have higher rates of returns and be more volatile. In a previous 
paper [24] the author gave a description of the joint dynamics of the market capitalizations in 
rank-based models, when the number of firms tends to infinity (see also ffTTTl for related results). 
Here, the corresponding limit is investigated in the context of volatility-stabilized models. 

The dynamics of the capitalizations in volatility- stabilized models is given by the unique 
weak solution to the following system of stochastic differential equations: 



dX.it) = ^S(t)dt + y/Xi(t)S(t) dWiit), 1 < i < N, (1.1) 

which is endowed with an initial distribution of the vector (Xi(0), . . . ,Xn(0)) on [0, oo)^. 
Hereby, rj is a real number greater than 1, S(t) = Xi(t) + ■ ■ ■ + X N (t) and W\, . . . , Wn is 
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a collection of N independent standard Brownian motions. We refer to section 12 of Q for a 
construction of a weak solution to (11.11 ) and an explanation why it is unique. 

We will analyze the limit of the path of empirical measures J2iLi ^Xi(t) corresponding to 
(11.11) . after a suitable rescaling of the time parameter, when N tends to infinity. The slowdown of 
the time by a factor of N is needed to observe a non-degenerate limiting behavior. Heuristically, 

this can be inferred from the appearance of the process S(t) — Xi(t) + h X N (t), an order 

N object, in the drift and diffusion coefficients of the processes X l5 . . . , X N . 

We show that the limit of the sequence of laws of $2 f=1 $Xi(t/N),N £ N exists and that un- 
der the limiting measure the degenerate linear parabolic equation (11.51) below must be satisfied 
in the weak sense with probability 1. Using a stochastic representation of the solution to (|1.5I ). 
we can determine the latter explicitly. Hence, our results allow to approximate the evolution of 
the capitalizations in a large volatility-stabilized market by the solution of the limiting equation 
(11.51) . Moreover, in the context of stochastic portfolion theory (see e.g. O, Q) one is inter- 
ested in the behavior of the rank statistics of the vector (Xi(t), . . . , X N (t)) of capitalizations. 
Since these are given by the jj, j^, . . . , ^-quantiles of the empirical measure jj Yn=i 5x t {t), our 
results can be also used to approximate the evolution of any finite number of ranked capitaliza- 
tions by the evolution of the corresponding quantiles of the solution to the partial differential 
equation (11.51) . In addition, the stochastic representation mentioned above shows that the so- 
lution to the equation (11.51) is given by the one-dimensional distributions of a time-changed 
squared Bessel process and, thus, establishes a new connection between volatility-stabilized 
market models and squared Bessel processes (see [d for further connections). The latter were 
analyzed in much detail in the works fl2TTl . [1221 and Il23l among others. 



Independently from the field of stochastic portfolio theory, systems of interacting diffusion 
processes play a major role in statistical physics. In particular, systems of diffusions interacting 
through their empirical measure (mean field) were studied in the literature by many authors, see 
e.g. HOl, ffHIl, [I9]], H3]], dU], HTJ. We remark that the system (O) can be cast into 
the framework of iTTOl . since the drift and the diffusion coefficients in the i-th equation of the 
system (11.11) can be expressed as functions of the empirical measure of the particle system and 
the position of the i-th particle. However, the generator of the particle system is not uniformly 
elliptic on [0, oo)^ and the same is true on [0, oo) for the elliptic differential operator on the 
right-hand side of the equation (I1.51 ). For this reason, the results of IfTOll do not carry over 
directly to our setting. Nonetheless, we adapt some of the techniques developed there to our 
case. 



The time- varying mass partition 

aAt) = , x Xi ® — , 1 < i < N (1.2) 

w x 1 (t) + --.+x N (ty - - 

is referred to as the collection of market weights in the mathematical finance literature and 
describes the capitalizations of the firms as fractions of the total capitalization of the market. 
The collection of market weights has the remarkable property of being extremely stable over 
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time for all major financial markets (see [|6) for plots of the market weights in multiple real- 
world markets). This was explained to a large extent in the context of rank-based market models 
in US and ||20t 

The model (11.11) incorporates the empirically observed fact that the capitalizations of firms 
with a small market weight tend to have a higher rate of growth and to fluctuate more wildly. In- 
deed, this becomes apparent from the dynamics of the logarithmic capitalizations corresponding 
to (PI): 

d(lo S Xi(t)) = dt + -;L= dWi{t), l<i<N (1.3) 

and the assumption 77 > 1. A detailed analysis of the evolution of the market weights under the 
model (11.11) and the corresponding invariant measure can be found in |fT9l . 

We assume the following condition on the initial values Xl(0), . . . , Xn(0) of the capital- 
izations. 

Assumption 1.1 The laws of — J2f=i ^(0). N G N on Mi([0, 00)), the space of probability 
measures on [0, 00) endowed with the topology of weak convergence of measures, converge 
weakly to S\for some X G Mi([0, 00)) with a finite first moment, the quantities E[5(0)] and 
K[S(0) 2 ] are finite for all N G N, and it holds 



JTF — m A, llm 



lim — — — = m\, lim — — — = m x , (1-4) 



where m\ = Jj ^ x X(dx). In addition, we make the non-degeneracy assumption m\ > 0. 

We remark at this point that Assumption II .H is, in particular, satisfied if the random vari- 
ables Xl(0), . . . , Xn(0) are i.i.d. and distributed according to a measure A with finite two first 
moments and m\ > 0. This is a consequence of Varadarajan's Theorem in the form of Theorem 
11.4.1 in 

In order to formulate our main results we introduce the following set of notations. We 
write Mi(R) and M 1 ([0, 00)) for the spaces of probability measures on the real line and the 
non-negative half-line, respectively. We metrize both spaces in a way compatible with the 
topology of weak convergence of measures. Moreover, for a positive real number T, we 
let C([0,T],Mi(R)) and C([0, T], Mi([0, 00))) be the spaces of continuous functions from 
[0,T] to Mi(R) and from [0,T] to Mi([0, 00)), respectively, endowed with the topology of 
uniform convergence. In addition, we introduce the time-changed capitalization processes 
Yi(t) = Xi(t/N), t > 0, 1 < i < N and the corresponding path of empirical measures 
Q N (t) = ■h Y^i=i $Yi(t), t G [0, T] on an arbitrary finite time interval [0, T], which is considered 
to be fixed from now on. Finally, we let be the distribution of the random variable g N (t), 
t G [0, T] on C([0, T), Mi( 



The main results of this paper are summarized in the following theorem. 
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Theorem 1.2 Under Assumption li.il the following statements are true. 

(a) The sequence Q?, N G N converges weakly to a limit Q™. Moreover, Qj? is a Dirac delta 
measure, whose unique atom q is given by the unique distributional solution of the Cauchy 
problem 

dg ri vt dg 1 m d 2 (xg) „ ^ x 

¥ = -2 e2mA ^ + 2 e2mA ^' (L5) 
g(0) = A (1.6) 

in C([0,T),Mi([0,oo))), where m x = J [0oo) x X(dx). 

(b) Let Z(t), t > be a squared Bessel process satisfying the stochastic initial value problem 

dZ(t) = 7}dt+ s/Z(t) d(3(t), t>0, (1.7) 
C(Z(Q)) = X, (1.8) 

where (3 is a standard Brownian motion and C(Z(0)) denotes the law of Z(0). Then the 
unique distributional solution to the Cauchy problem (II .51) . (11.61) inC([0, T], Mi([0, oo))) is 

given by the one-dimensional distributions of the time-changed process Z y J * e ris l 2 m\ dsj, 

t G [0, T\. 

(c) Let g(t), t G [0,T] be the only atom of the measure Q^- Then for every t G (0, T] the 
measure g(t) is absolutely continuous with respect to the Lebesgue measure Leb on [0, oo) 
and the corresponing density is given by 

where J(t) = f Q e^^mx ds and is the Bessel function of the first kind of index r\ — 1. 



Remarks. 

( 1 ) By a distributional solution of the problem (11.51) . (|1.61 ) in C( [0, T] , Mi ( [0, oo) ) ) we mean an 
element g of C([0, T], M x ([0, oo))) which satisfies the system (ITT71) . (ITT81) . where S(R) 
is the space of Schwartz functions on R and (7, /) denotes J R f dj for any / G «S(R), 
7 e Mi(M). 

(2) In E4l the limiting dynamics was derived (under some assumptions) for a different class of 
interacting diffusion processes, which go by the name of rank-based models in the context 
of stochastic portfolio theory. There, the limiting equation for the cumulative distribution 
function of the empirical measure of the logarithmic capitalizations was given by the porous 
medium equation, that is, a nonlinear non-degenerate parabolic partial differential equation. 
It was also shown there that its weak solution w can be represented by the one-dimensional 
distributions of the process with the dynamics 

dX(t) = /j(w(t,X(t))) dt + a(w(t,X(t))) d/3(t), (1.9) 
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where [i and a are functions depending on the parameters of the model. In contrast, the par- 
tial differential equation (11.51) is a linear degenerate parabolic differential equation, which 
admits the stochastic representation of Theorem II .21 (b). Hence, although the rank-based 
and the volatility- stabilized models share multiple common features (such as the monotone 
dependence of the drift and diffusion coefficients of the logarithmic capitalizations on the 
market weights), their limiting behavior differs significantly. Indeed, it was shown in HI 
that the weak solution of the porous medium equation may fail to be differentiable in the 
spatial variable for all t > in contrast to the findings in Theorem ll .21 (c) for the equation 
(11.51) . In addition, a big difference from the applicational point of view is the explicitness 
of the solution to the equation (11.51) . whereas (in general) no explicit formula for the weak 
solution of the porous medium equation is known. It is also remarkable that with the usual 
parametrizations of the two models as in [7 |, the time in the volatility- stabilized models 
(11.11) has to be slowed down by a factor of N to observe non-degenerate limiting behavior, 
whereas this is not the case in the rank-based market models. 

(3) Let d be a metric on Mi (R) which metrizes the topology of weak convergence of probability 
measures (such as the Levy metric or any other metric in section 1 1.3 of [5]), and let g?[o,t] 
be the metric on C([0, T), Mi(R)) given by 

d[o,n(&>6)= SU P &(*))■ (i-io) 

te[o,T] 

Then dmr] makes the space C([0, T], Mi(R)) into a separable metric space (see e.g. Theo- 
rem 2.4.3 in [25 1). Moreover, combining part (a) of Theorem II .21 with problem 6 in chapter 
9.3 of [5], we conclude that the sequence g N , N E N converges to the path of measures g 
of Theorem 0(c) in probability on (C([0, T], Mi(R)), d [0 , T ]), that is, 

Ve > : lim P( sup d(g N (t), g(t)) > e) = 0. (1.11) 

N ^°° te[o,T] 

This gives an alternative way of stating part (a) of Theorem II .21 

(4) The transition densities of the squared Bessel process Z in Theorem 11.21 (b) are known 
(see e.g. Corollary 1.4 in chapter XI of [23 1), so that part (c) of Theorem 11.21 is a direct 
consequence of parts (a) and (b) of Theorem 11.21 

The rest of the paper is organized as follows. Part (a) of Theorem II .21 is proven in sections 
2.1 and 2.2. Its proof is divided into three parts. Firstly, in Proposion l2.1l it is shown that the 
sequence Qj,, N 6 N is tight. Its proof relies on the characterization of compact subsets of 
C([0,T],Mi(R)) obtained in Qil and a characterization of compact subsets of C([0,T],R), 
the space of continuous real- valued functions on [0, T] endowed with the topology of uniform 
convergence, given in [26J . Secondly, in Proposition l2.3l we prove that under every limit point of 
the sequence , N G N the Cauchy problem (11.51 ). (11.61 ) is satisfied in the distributional sense 
with probability 1. Here, we use arguments from the theory of convergence of semimartingales 
in the spirit of lfT2l . The main challenge in these two parts is to deal with the unboundedness 
of the drift and diffusion coefficients in the dynamics of the processes Yi, . . . , Yjy. Thirdly, in 
Proposition l2.4l we demonstrate that the problem (11.51 ). (|1.6I) has a unique distributional solution 
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in the space C([0, T], Mi([0, oo))). This is achieved by transforming the uniqueness problem 
into an existence problem and by applying existence results on boundary value problems for 
non-degenerate linear parabolic equations on bounded subsets of [0, T] x R, on which the dif- 
ferential operator of equation (11.51 ) is uniformly parabolic. After that, we give the proof of part 
(b) of Theorem II .21 in section 2.3 using methods of stochastic calculus. 



2 Law of large numbers 
2.1 Tightness 

In this subsection we will combine the characterization of relatively compact sets in the space 
C([0, T), Mi(R)) of 0H with a characterization of relatively compact subsets of C([0, T), R) 
in Il26l to prove the tightness of the sequence Q?, N G N. 



Proposition 2.1 The sequence Q%, N eNis tight on C([0, T], M 

Proof. 1) Let C C (R) be the space of compactly supported continuous functions on R endowed 
with the topology of uniform convergence. We fix an arbitrary e > and a countable dense 
subset {/i, / 2 , . . . } of C C (R) such that each f r is twice continuously differentiable. Moreover, 
for every 7 G Mi(R) and every function / on R, which is integrable with respect to 7, we write 

(7, /) for / R / d'y. 

From the proof of Lemma 1.3 in [fTOl we see that it is enough to find a compact set K in 
Mi(R) and compact sets K u K 2 , . . . in C([0, T], R) such that for all iV G N: 

Q^({£ G C([0,T],M 1 (R))|Vt G [0,T] : £(t) G K }) > 1 - e, (2.1) 
Q^({e G C([0,nMi(R))|(£(.),/r) G ^ 1 - ^2" r , r > 1. (2.2) 

To define -Ko we introduce the function tp(x) = \x\ and use the non-negativity of the pro- 
cesses Yi, . . . , Y N together with the dynamics of the processes X\,..., X N to conclude 



i=l 



where B,(t) = N l l 2 Wi(t/N), 1 < i < N, S Y (t) = Y 1 (t) + - ■ -+Y N (t). From the representation 
of Xi, . . . ,Xjy as time-changed squared Bessel processes (see equations (12.7)-(12.9) in Q) 
and remark (ii) after Corollary 1 .4 in chapter XI of [1231 (note that their dimension parameter 5 
corresponds to our 2rf) it follows that the return time to of the processes Xi, . . . , X N is infinity 



rt y/Yj(s) 

and a standard Brownian motion by Levy's Theorem. As a consequence we have 



with probability 1. Hence, the process B(t) = J2i=i Jo 1 ^ dBj(s), t > is well-defined 

v S (s) 



d(g N (t), <p) = % N (t), <p)dt + { -^%^-dB(t). (2.3) 
^ V A 
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The latter equation is a Black-Scholes stochastic differential equation and it is well-known that 
its unique strong solution is given by 

(Q N (t), <p) = (Q N (0), <p) exp {{rj/2 - (2W)" 1 )* + N^Bit)) . (2.4) 

Thus, for every C > and all iV G N we have 

>p (g N (t),<p))>C 

0<t<T 

< pf(^(0),^)exp(A^- 1/2 sup B(t)) > Cexp(-(r]/2- (2N)- 1 )T) 

\ 0<t<T 



A routine computation involving Chebyshev's inequality and (11.41) shows that the sequence of 
random variables (^(0), <p) exp(iV~ 1//2 sup 0<t<T B(t)), N G N converges in probability to the 
constant rn\. Hence, by choosing C large enough, one can make the latter upper bound smaller 
than e for all N G N. Thus, we can let K be the closure of the set 

{7GM 1 (M)|(7,^)<C} 

in Mi(R), which is compact by Prokhorov's Theorem. 

2) To prove the existence of the sets K x , K 2 , . . . with the desired properties it suffices to show 
that for any fixed r G N the sequence of probability measures Q" Jr , N G NonC([0,T],l) 
induced by Q?, N G N through the mapping £ !->■ (£(.), f r ) is tight. To this end, we fix an 
r G N and aim to deduce the tightness of the sequence Q^'^ , N E N from Theorem 1.3.2 of 
Il26ll . To do this, we need to show 



lim inf Qr' Ir (\y(0)\ < 



and 



VA > : limlimsupQ^' /r ( sup \y(t) - y(s)\ > A) = 0, 



C40 jv-kx> V ( s ,t)eA c 

where = {(s,t)\0 < s < t < T,t — s < (}. The first assertion follows immediately by 
considering 9 > sup,. eK \f r (x)\. 

To prove the second assertion, we first rewrite it in terms of X\ 1 . . . , X N : 



1 N 

VA > : limlimsupPf sup - V(/ r (X,(t)) - f r {Xi(s))) 
C4-0 Af^oo V ( S! t)eA (tN iV 1 r— : 



N 

' :> A 

i=l 



where v4 c>A r = {(s,t)|0 < s < t < T/N, t - s < C/N}. Next, we apply Ito's formula to the 
process jj Yli=i fr{Xi{t)), t > and conclude that it holds 



1 N 

- fr{Xi(t)) = D(t) + M(t), t > 0, (2.5) 



i=l 



where 

N 



D ^ = ^ E / lf'r(Uu))S{u) + ^'(X^X^Siu) du, (2.6) 

l N r* 

M W = JJzZ K(Xi(u)WXi(u)S(u) dWi(u). (2.7) 
1=1 ^° 

Thus, for every fixed A > the latter limit can be bounded from above by 

limlimsupPf sup \D(t) - D(s)\ > A/2) 

CIO at^oo \ o<s<t<T/N,t-s<(/N ' 

+ limlimsupPf sup \M(t) - M(s)\ > A/2) , 

C4-0 n^oo \ 0<s<t<T/N,t-s<(/N ' 

which we will call expression (*). We will show that the first summand in (*) is zero in step 3 
and that the second summand is equal to zero in step 4. 

3) To bound the first summand from above, we set 

R = max(77/2sup \f r (x)\, l/2sup \f?(x)\) (2.8) 

xeR xeR 

and use the definition of the process D to make the estimates 

Pf sup \D(t) -D(s)\ > A/2) 

< P( sup / S(u) + S{u) 2 /N du > A) 

^ 0<s<t <T/N,t~s<(/N Js 2it / 

< P( sup (S(t) + S(tY/N) > £t) 

V 0<t<T/N ^-KC, 7 

= P ( sup S (t)>- 1 + ^ + 2A <W N ). 

^0<t<T/N 2 / 

To estimate the latter upper bound further, we use the dynamics of the processes X 1: . . . , X N 
to find 

dS(t) = S(t)^-dt + V x i(t)S(t)dWi(t) = S(t)^-dt + S(t)dB(t), 

1=1 

where B{t) = ^2f =1 J* * dWi(s), t > is a standard Brownian motion due to the same 

argument as for the process B in step 1 of this proof. Thus, S satisfies the Black-Scholes 
stochastic differential equation and is given explicitly by 

S{t) = S(0) exp ((77JV/2 - 1/2)* + ,t>0. (2.9) 
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Hence, the latter upper bound is not greater than 

p(f)exp((, J V/2-l/2)r/ J V + 8U p g M )> ±« ), 

V iV V Q<t<T/N ' I J 

From (11.41) and Chebyshev's inequality it follows that the sequence of random variables 

S(0) 



■ exp ((r]N/2 - \/2)T/N + sup B(t)) , N G N 

A ^ 0<t<T/N ' 

converges to the constant mxe nT ^ 2 in probability in the limit N — > oo. Thus, the latter proba- 
bility converges to in the limit iV — > oo for all ( small enough. 

4) To show that the second summand in expression (*) is zero, we first note that for every pair 
< s < t < T/N with t - s < C/N there is a k G N with s, t G [k(/N, (k + 2)(/N]. We use 
this observation and the union bound to conclude 

P( sup \M(t) -M(s)\ > A/2' 

^ 0<s<t<T/N,t-s<C/N ' 
LT/CJ-1 

< V PI sup |M(t)-M(fcC/iV)| > A/4). 

V kC/N<t<(k+2)C/N ' 

where [ J denotes the function taking a real number to its integer part. 



(2.10) 



Next, we use (12.91) to compute 



E[S{t)\ = E[S(0)]exp(rjNt/2), 
E[S(t) 2 } = E[S(0) 2 ] exp((r/iV + l)t). 



The inequality /;(X i (t)) 2 X i (t)S(t) < su Pa:6K \ft(x)\ 2 S(t) 2 , Fubini's Theorem and d2~T0l im- 
ply that the process M(t), t > is a martingale. Applying the L 2 -version of Doob's maximal 
inequality for non-negative submartingales we obtain 



IT/CJ-1 



V p( sup |M(t) -M(fcC/iV)| > A/4) 



fc=0 >K/N<t<(k+2)C/N 

LT/CJ-i 



< & E E[(M((fc + 2)C/iV)-M(K/iV)) 2 ]. 
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fc=0 

By the Ito isometry the latter expression can be computed to 
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LT/CJ-l iV r (k+2)C/N 

V VE fXXi(u)fXi(u)S(u) du 

t~n ~1 l JkC/N 



A 2 N 2 

k=0 i=l 



32 SUp eK |/r(-'c)|' „ / , 2 . 

" A^iV 2 " ' ' ( ' du 



(IT/Q+1K/N 



By Fubini's Theorem (note that the integrand is non-negative) and (12.101) one deduces that the 
right-hand side is equal to 



32su PxgR |/;(x)| 2 E[g(0)- 
A 2 N 2 (r]N + 1) 



e W ((r,N + l)([T/C\+l)C/N 



Using (11.41) we conclude that the latter expression tends to in the limit N — > oo for any fixed 

C>o. □ 



2.2 Identification of the limit point 

In this section we will uniquely characterize the limit points of the sequence , N G N and 
thereby complete the proof of Theorem II .21 (a). To this end, we fix a convergent subsequence 
Q^ k , k G N of the sequence Q% , N G N and let be its limit. By the Skorokhod Representa- 
tion Theorem in the form of Theorem 3.5.1 in 0], there exist C([0, T], Mi (R)) -valued random 
variables Qk, k G N with laws Q^ k , k G N converging to a limiting random variable with 
law almost surely. We start with the following lemma. 

Lemma 2.2 The sequence of functions t \-t (Q k (t), x), k G N converges in the space C([0, T], K.) 
tot ^ m\e r,t l 2 in probability in the limit k — > oo. 



Proof. First, we fix a G N and with a minor abuse of notation write S(t) for Xi(t) + ■ ■ ■ + 
X Nk (t). By equation ( 12.91) , we have for every e > 0: 



pf sup (g k (t),x) -m x e^ 2 

\ 0<t<T 



> £ 



Pf sup exp (jj/2 - l/(2JV*))t + iV; X/ V(t)) - m A e^ /2 



> e 



S(0) 



< Pf sup exp f - */(2JV fc ) + A^V^)) - m A > ^e^ 2 ) 



1/2 — 

where W(t) = B(t/N k ), t > is a standard Brownian motion. By Girsanov's Theorem 

the process exp ( — t/ (2N k ) + N^ 2 W(t) J — m x , t > is a martingale. Hence, we can 

apply the L 2 -version of Doob's maximal inequality for non-negative submartingales to estimate 
the latter upper bound from above by 



-E 



exp ( - T/(2N k ) + N~ 1/2 W(T)) 



2n 



,rfT 



-E 



E[S(0)1 2 



nsm eT/Nk _ 2mx nsm +ml 



N, 



It follows from (11.41) that the latter expression tends to in the limit k — V oo. This finishes the 
proof of the lemma. □ 



10 



Now, we are ready to prove that the system (11.51 ). (11.61 must hold in the distributional sense 
almost surely under any limit point of the sequence , N G N. 

Proposition 2.3 Let be the limit of a convergent subsequence Q^ k , k G N of the se- 
quence Qt, N G N. Then under the system (11.51) . (11.61) is satisfied in the distribu- 
tional sense almost surely. Moreover, if is a random variable with law Q™, then it holds 
Qoo G C([0, T], Mi([0, oo)) with probability 1. 

Proof. 1) Fix a as in the statement of the proposition. Let {g 1 , g 2 , . . . } be a dense sub- 
set of the space <S(R) of Schwartz functions on R with respect to the topology of uniform 
convergence of functions and their first and second derivatives, such that each g r is infinitely 
differentiable with compact support. We claim that in order to prove the first assertion of the 
proposition it suffices to show that 



holds for all r G N and t G [0, T] with probability 1 under . Indeed, this would imply that 
the system (12.111) . (12.121) is satisfied for all g G «S(R) and all t G [0, T] with probability 1 under 
Q^. This would yield the first assertion of the proposition. 

2) Since a countable union of null sets is a null set, it is enough to show that the system (12.1 II) . 
(12.121) is satisfied for a fixed function g r and all t G [0, T] with probability 1 under Q™- It is 
clear that equation (12.121) is satisfied with probability 1 due to Assumption ll.li In order to show 
that equation (|2.11l) holds, we use Ito's formula to compute 




(2.11) 



g(0) = A 



(2.12) 



i=l i=l 



i=i 



Thus, one has the dynamics 



d(g (t),g r ) = D r (t) dt + 



^ m Y,g' r {Y l {t))^it)S Y it)dB l {t), (2.13) 



i=i 



where 



i=l 




11 



t > 0, and Bi(t) = N l / 2 Wi(t/N), 1 < i < N, S Y (t) = Y 1 (t) + • ■ ■ + Y N (t) as before. 
Moreover, the inequality 



/ r (Ximy/^WW) < snp\g' r (x)\S Y (t), 



(2.14) 



xGM. 



Fubini's Theorem and (12.101 ) show that the process 



M r( f ) = / 9r{yi{n)Wy i {u)s Y {u)dB i {u), t > o (2.15) 

1=1 J° 

is a martingale. Hence, using the L 2 -version of Doob's maximal inequality for non-negative 
submartingales and the Ito isometry, we obtain for every e > 0: 



F( sup (g N (t),g r )-(g N (0),g r )- [ D r {u) du 



• 0<t<T 

< e- 2 E[M r (T) 2 } 

N rT 



> e 



e 2 N 3 



i=i 



J> / g' r (Y i (u)) 2 Y i (u)S Y (u)du 



< 



sup 



Wr(x)\ 2 



e 2 N 3 



E 



S y (w) 2 rfw 



Using S Y (u) = S(u/N), u > 0, Fubini's Theorem and (12.101) we can compute the latter upper 
bound to 

su Px . eM K(x)| 2 iVE[S(0) 21 



e 2 N 3 



T)N + 



HA (^,(r,N+l)T/N _ ^ _ 



(2.16) 



This expression tends to in the limit iV — > oo for every e > due to (11.41) . 



3) Next, we recall the definition of the random variables g k , k £ N and g^ prior to Lemma I2i2l 
In view of the latter, we may and will assume that the sequence of functions t !->■ (gk{t), x), 
k G N converges to t mxe^^ 2 in the space C([0, T],R) with probability 1 (otherwise we 
pass to a suitable subsequence). It follows that the random variables E k given by 



sup 

0<t<T 



[Qk(t),9r) ~ (Qk(0),9r 



[g k (u),r]/2g! r + x/2g") ■ (g k (u),x) du 



converge almost surely in the limit k — >■ oo to 



sup 

0<t<T 



(t),9r) ~ (0oo(O),</ r ) 



which we call Soo. Finally, using the Portmanteau Theorem and the final result of step 2 we 
obtain for every e > 0: 



P(S co > e) < liminf P(S fc > e 



k— >oo 



lim inf P I sup 

fc->oo V <f<T 
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DJu) du 



> e 



0. 



Since the law of is given by , it follows that equation (12.111) holds Q~ -almost surely. 

4) To prove the second assertion of the proposition, we note that £>fc(£)([0, 00)) = 1 holds 
for all k G N and t G [0, T] almost surely. This is a consequence of the representation of 
Xi, . . . , Xn as time-changed squared Bessel processes (see equations (12.7)-(12.9) in [7J) and 
the properties of the latter (see e.g. chapter XI of [23]). Thus, the Portmanteau Theorem implies 
f?oo(*)([0, 00)) = 1 for all t G [0, T] on the same set of full probability. □ 

In view of Propositions l2Tl and l2.3[ the proof of Theorem 1 1 . 2 1 (a) is complete, once we show 
that the solution of the Cauchy problem 

VgeS(R): (g(t),g) - (g(0), g) = m x j (e(s),|</ + ~ V) ds, (2.17) 
q(0) = A (2.18) 
in C([0, T], Mi([0, 00))) is unique. 

Proposition 2.4 The solution of the Cauchy problem (12.171) . (12.181) in the space C([0, T] , M x ([0, oc 
is unique. 

Proof. 1) Let /i, v G C([0,T], Mi([0, 00))) be two solutions of the problem (12T7T) . (l2T8i 
Moreover, define the operator 



dh rjtridh ntxd 2 h 

— m\e 2 — — m\e 2 — -— : 

dt 2dx 2dx 2 



yy lb 'I'll i_/ 1 u tXj kJ I L 

Lh = — rn x e^-- m x e^ — — — - (2.19) 



acting on the space C^ ,2 ([0, T] x R) of continuous real-valued functions on [0, T] x R having 
compact support in (0, T) x R, a continuous time derivative and two continuous spatial deriva- 
tives. For the same reason as in the remark preceeding Theorem A. 1 in IfTOll . the problem (12. 171 ). 
(12.181) is equivalent to the problem 

\/h G C?([0, T] x R) : (g(t), h(t, .)) - (f)(0), h(0, .)) = f\ 6 (s), (Lh)(s, .))ds (2.20) 

Jo 

with initial condition g(0) = A. This follows by an approximation of functions h G C]' 2 ([0, T] x 
R) by functions on [0, T] x R which are Schwartz functions in x for every fixed t and piecewise 
constant in t. Since for all h G C c 1>2 ([0, T\ x R) it holds (/i(T), h(T, .)) = (i/(T), h(T, .)) = 0, 
we conclude from (12.201 ): 



pT pT 

y1,2/ 



Vh G C c ll2 ([0,T] x R) : / (fi(s), (Lh)(s, .)) ds = (u(s), (Lh)(s, .)) ds. (2.21) 

Jo Vo 

2) In view of the latter equation, it suffices to show that the space 

Lh, h G C c 1,2 ([0,T] x R) 

contains all functions on [0, T] x R of the form (t, x) 1— >■ l(t)r(x), where Z is an infinitely 
differentiable function on [0, T] with compact support contained in [to, t%\ for some < to < 
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ti < 7 and r is an infinitely differentiable function on K. with compact support contained in 
[a, b] for some < a < b. To this end, we first let 72 be a bounded closed rectangle of the form 
[to, ti] x [a, 6] which is contained in (0, T) xl and such that [t , ti] x [a, 6] is contained in the 
interior of 72.. By Theorem 9 in section 1.5 of [8] there exists a function h : 1Z — > R such that 
the derivatives ||, |4| exist and are continuous on 72 and it holds (Lh)(t, x) = l(t)r(x) in 
the strong sense on 72. 

Next, we apply the maximum principle in the form of Lemma 4 in section 2. 1 of [H on the 
rectangles 

[to,?i] x [a,a], x [a, 6], fF ,^i] x [ML fo,*o] x [a, 6] 

to the functions h and — h to conclude that the function h is constant on the complement of 
[to, ti] x [a, 6] in 72. Hence, without loss of generality we may assume that h is constantly equal 
to zero on the complement of [to, ti] x [a, b] in 72 (otherwise we subtract a constant from h). 
Thus, we can extend h to a function h 6 C^' 2 ([0, T] x K) such that h = h on 72 and h = 
on the complement of 72. It follows that (Lh)(t,x) = l(t)r(x) holds in the strong sense on 
[0, T] x R, which finishes the proof. □ 

2.3 Stochastic representation of the limit point 

In this subsection we prove part (b) of Theorem [L2l 

Proof of Theorem\L2\(b). In view of Proposition [23] it suffices to show that the one-dimensional 
distributions f(t), t G [0,7] of the process T(t) = z( J^e^mx ds\t e [0,7] form a dis- 
tributional solution to the Cauchy problem (fT3T> . (fL"6T ) in C([0, T], Mi([0, oo))). It is clear 
that t G [0, T] satisfies (fL6l) and that it is an element of C([0, T], Mi([0, oo))), since the 
squared Bessel process Z has continuous paths and takes values in [0, oo) (see e.g. chapter XI 
in [|23]|). To prove that £(t), t G [0, 7} satisfies (11.51) . we use the time-change formalism for 
Brownian motion (see e.g. lfT8l . chapter 8.5) to deduce 



= e nt/2 m x ^ dt + e^ /4 v^U\/f(t) dp{t), 
t > 0, where is an appropriate standard Brownian motion. Hence, Ito's formula shows 

g(T(t)) - g(T(0)) = f ™ x e"l* (|</(r(s)) + l -Y{s)g"{V{s))) ds 

+ f ' g\T(s))e^y/^y/f(s)dP( S ) 
Jo 

for all g G S(R) and t G [0,7). Moreover, the expectations E[r(t)], t G [0,7] are finite and 
uniformly bounded (this is evident from the definition of squared Bessel processes in chapter 
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XI of Il23ll if 2r] is an integer; in the general case, this can be deduced by comparing Z with a 
squared Bessel process of integer index 2r)' > 2rj using, for example, the comparison theorems 
of section 3 in chapter IX of ll23~T0 . Thus, by taking the expectation in the latter equation and 
applying Fubini's Theorem we get 

(m,9) - (£(0),<?) = / (as),m x e^ 2 (^g'+ -xg")) ds (2.22) 

for all g E S(R) and t E [0, T). Thus, £{t), t E [0, T] solves the equation (O) in the distribu- 
tional sense. □ 
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